Narrow gap Luttinger liquid in Carbon nanotubes 
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Electron interactions reinforce minigaps induced in metallic nanotubes by an external field and 
turn the gap field dependence into a universal power law. An exactly solvable Gross-Neveau model 
with an SU{A) symmetry is derived for neutral excitations near half-filling. Charge excitations, 
described by a sin-Gordon perturbation of Luttinger liquid theory, are composite solitons formed 
by the charged and neutral fields with two separate length scales. Charge compressibility at finite 
density, evaluated in terms of inter-soliton interaction, exhibits a crossover from overlapping to 
non-overlapping soliton state. Implications for the Coulomb blockade measurements are discussed. 



Electron interactions create a peculiar strongly corre- 
lated ID electron system in metallic Gabon nan- 
otubes, the thinnest and the cleanest among the currently 
available nanoscale quantum wires. Luttinger liquid the- 
ory of nanotubes predicts that, since tube diameter 
is larger than Carbon separation, the ID electron cou- 
pling is mainly accounted for by the long-range electron 
interaction (forward scattering), while the exchange and 
Umpklapp scattering, as well as backscattering, are rela- 
tively weak Recent experimental work focused 
on Luttinger liquid effects in tunneling, observed as char- 
acteristic power laws in the tunneling current dependence 
on bias voltage and temperature. 

Here we discuss Luttinger liquid effects in nanotubes 
with a niinigap at the band center induced by an ex- 
ternal perturbation. Such a minigap can be opened by 
parallel magnetic field or by the intrinsic curvature 
of the tube . The field- and curvature-induced gaps 
were observed experimentally [ pljT^ ] and found to be in 
agreement with the noninteracting electron model. We 
show that electron interaction enhances the charging gap 
and makes it a power law function of the bare gap. 

This provides a unique situation, not available in other 
quantum wires, when Luttinger liquid effects are mani- 
fest in thermodynamical properties. For instance, in a 
gapped state induced by magnetic field, the bare gap is 
determined without any fitting parameters by Aharonov- 
Bohm flux through tube cross section, while the charg- 
ing gap is directly measurable via Coulomb blockade, a 
prominent feature of transport in nanotubes |l^-p^ . In 
the strong forward scattering limit, the power law rela- 
tion of the charging gap and magnetic field is character- 
ized by a universal exponent 4/5. We emphasize that 
the charging gap measurement is qualitatively different 
from the tunneling current measurement because it can 
be performed in thermodynamic equilibrium. 

Also, in the presence of interactions the charging gap is 
much enhanced compared to the neutral excitation gap, 
while in a noninteracting system the two gaps are pre- 
cisely equal. Large energy separation of the charged 
and neutral sectors results in high symmetry of the 
states in the neutral sector described by multiplets of 



the group 0(6) ~ 5'C/(4) derived from exactly solvable 
Gross-Neveau model. This picture, demonstrated in the 
situation when the gapped state is created by external 
field, is realistic, since intrinsic interaction-induced gap 
in nanotubes are believed to be extremely small 

Electron bands of metallic tubes form two pairs of 
spin-degenerate right and left branches intersecting at 
the band center [ pO[ . The Hamiltonian in the forward 
scattering approximation ||l|,|^ has the form 

where tpjix) is a two component wavefunction, p{x) = 

Sj=i V'/(3^)^i (^) is charge density, and v is Fermi ve- 
locity. The form of the forward scattering amplitude 
in Eq.(|l|) depends on the electrostatic environment. In 
a nanotube of radius r, in the absence of screening, 
V{q) = JV{x)eJi''dx = e^\n[{qr)-^ + l]. The substrate 
dielectric constant e reduces V{q) by a factor 2/{e+ 1). 

Field-induced gapped state ||Jl^ is described by 
adding to the Hamiltonian (|l|) a backscattering term 

Vcxt ^AoJ E ■=i^>iV'jrfa; . (2) 

In the absence of interactions, V{q) ~ 0, electron spec- 
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trum is e{p) — ± {y'^p^ + ^0) with the value of Aq de- 
pending on the backscattering mechanism. The magnetic 
field-induced gap H is linear in the field: Aq = hvcjy/r, 
where (f> — Trr^B/^o is the fiux through tube cross- 
section scaled by $0 = hc/e. For typical tube radius 
r ~ 0.5 nm and B ~ lOTesla, the gap Aq ~ lOnieV. 

We bosonize the Hamiltonian H = Ho + Vcxt in the 
standard way, using ipj cx e*^*^/r^/'^. The Gaussian 
part Tio is diagonalized using linear combinations of the 
bosonic fields = {e'j^'a) with 

eo = i(l, 1,1,1), ei = 1(1,-1,1,-1), 
62 = 5(1,-1,-1,1), 63 = i(l, 1,-1,-1) (3) 
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In this notation the Gaussian part of the Lagrangian de- 
scribes one charged and three (neutral) flavor modes: 

+ \j da:ELi(^M</'a)', K,^l + ^V{q) (4) 



[Ti = V ^ 1). Bosonizing Vcxt = Aq / dxY^'^-^{R^Lj 



h.c), we have 



Lx = -2A / dx 



47r$, 



(5) 



with A — Ao/r. The total Lagrangian L = Lq + C\ dis- 
plays the fundamental U{1) x SU{A) symmetry playing 
a major role in our analysis. Forward scattering makes 
the perturbation even more relevant. At zero chemi- 
cal potential the coupling A grows under renormalization 
and opens spectral gaps. 

There are several possible regimes in which a nanotube, 
described by Eqs.(^,(||), may exist. First, there is an 
insulating regime with the density at half-filling, where 
all excitations are gapped. Second, there are conducting 
states which can be realized by applying various exter- 
nal fields. These fields may close some gaps or even all 
of them, provided their magnitudes exceed certain crit- 
ical values. For example, by varying chemical potential 
one can close all the gaps and make the perturbation A 
irrelevant. This will lead to a transition into a metal- 
lic (Tomonaga-Luttinger liquid) regime. Fields breaking 
the SU (4) symmetry will not affect the charge sector and 
thus leave the system insulating, but may close some of 
the gaps in the flavor sector. 

In a nanotube near half-filling, with no external fields 
present except A, the perturbation (^) grows under RG as 
= {V where a = R-^'"^. By the self- 
consistency argument A(/') = the charge and flavor 
gaps are related to the nonintcracting gap Aq as 

Afl~ Ao(I?/Ao)(i-")/(5-"), Ach =^ /^i/'Afl, (6) 

where D = hv/r is ID bandwidth. Since Aq is propor- 
tional to external magnetic field, Eq.(|6|) predicts a power 
law scaling of the gaps versus an experimentally control- 
lable parameter. For high charge stiffness Kq 3> 1, the 
gap scaling exponent is universal, with the value 4/5. 

The separation (||) of the charge and flavor sector en- 
ergy scales enables one to integrate out the fast mode (f)Q 
in adiabatic approximation. This will generate an effec- 
tive action for the flavor modes with SU(4) symmetry. 
Separating 00 in Lagrangian (|^) one obtains 



(7) 

(8) 



C\ — — 4A j dx yui cos (i>a + U2 sin 



where 0i = •^/TT'/'i ■ Treating the slow fields Ui as adiabatic 
parameters and shifting the variable 



+ 77 = TT ^/^tan ^(712/^1) (9) 
we transform the total Lag 



£ = £o[0o + 7?] - / da;M[(/)a] cos^o 



(10) 
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with M[0a] = 2A (^1 + E^_^^ cos20a cos20fcj From 

Eq.(|lO|) we derive an effective action for the flavor sec- 
tor valid for energies well below the charge gap (|^) by 
integrating over the fast mode 0o- Let us examine the 
results of this integration. Writing the first term in ( [l0| ) 
as /:o[0o] + -CoM -h E(^,,{'?^'o}g(w^ Kqq^){q\q we note 
that the ry-dependent terms are strongly irrelevant. For 
example, £o[7?] contains squares of gradients 



9^01 sin 202 sin 203 -t- permut. 

1 + Hai^h cos 20a cos 206 



(11) 



consisting of series of operators with the minimal scaling 
dimension 2. The relevant contribution arises from the 
last term of Eq.(|l^). Integrating it over 0o we obtain the 
ground state energy of the sine-Gordon model 

_^^2/(2-a/4) j^^j ^ ^Q^g^ _ y^Y^^^^ 20, COS 20b + ... 

= "S'EaT^fc ■ cos20acos20h : +... 

where 5 is a suitably renormalized coupling constant and 
dots stand for less relevant operators. The normal order- 
ing is taken with respect to the new cut-off G. 

Thus at energies smaller than the cut-off we obtain 
the effective action for the flavor modes in the form of 
the 0(6)-- SU(4) Gross-Neveau model: 



da; 



1, 



2Ea=l(^M'/'a) -5Ea#b ■ COS 20a COS 2 0h 

= / da; {iXili^^i^Xj ~9 ■ (XjXj)(XfeXfe) :] (13) 



where Xj U — 1j ...j6) are Majorana fermions. The lat- 
ter model is exactly solvable, the spectrum consists of 
3-1-6-1-3 relativistic particles with masses m, y/2m, m 
transforming according to different representations of 
the 0(6)~ SU(4) group Their dispersion law is 

Ei{p) = {v^p^ -\- ra^Y^"^. The Majorana fermions them- 
selves belong to the vector representation of the group 
and carry mass V2m. The mass gap m depends on the 
bare parameters and can be crudely estimated as follows. 
The cut-off G is of order of the charge gap Ach; it may 
be somewhat greater than A^h because the charge mode 
has a large velocity Vc = K^^^v ^ v and therefore can 
be considered as fast in a larger region of the momentum 
space than for Vc = v. The value of the coupling constant 
g is such that the mass of the Gross-Neveau model ( p^ ) 
coincides with the physical spin gap (O) : 
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G(Ao)e- 



-27r/3(; 



Afl = /(if)Ao(D/Ao 
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(14) 



At present we cannot estimate the function f{x). 

As for the charge excitations, they carry exactly the 
same quantum numbers as an electron. This happens 
because whenever a soliton of (j)o is created, it acts as an 
effective potential for flavor fields forming bound states. 
The flavor fields coupled to a (/)o-soliton give it the corre- 
sponding quantum numbers. Therefore, charge and fla- 
vor cannot be separated from each other in the charge 
sector. 

In the sin-Gordon problem electron is repre- 

sented by a soliton of one of the fields <i>j . Soliton spa- 
tial size can be estimated from a variational principle. 
Due to large charged stiffness K, the elastic energy of 
the soliton is dominated by the field (j)o. According to 
and (^), the soliton energy, estimated from the en- 
ergy in the space interval / where the field ^jix) varies, 
is E{1) ~ [Kqi^il^^ + '^O ^ with renormalized coupling 
X'{1) = X{r/l)'^^+°'y^. Soliton size w can be found by 
minimizing the energy with respect to I. Ignoring the log- 
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arithmic /-dependence of K, this gives w ~ {K/\'{w)) ' . 
The soliton energy E{ii}) ~ {K\' [w))^^"^ coincides with 
the charge gap Ach in (||). 

A peculiar feature of the composite charge soliton is 
the presence of two different length scales, because the 
flavor fields 0i,2,3 vary faster than the more stiff charge 
field (pQ. Let us consider a variational solution of the 
problem £o + -Ca with one of the fields $j varying be- 
tween two minima of the energy (^. To be specific, we 
consider a soliton of the field <I>i = ^ (00 + ••• + '/'s) in 
which $1 changes by ^/tt, while '&i,2,3 do not change. In 
this case all fields 0o,...,3 change by 7r/2. The length w 
over which the field c/io changes is much larger than that 
for (/)i,2,3- Thus the variational problem for 0o can be 
treated in a 7r/2 step approximation for 0i,2,3, as 



\K{<1)'qY = 47rA' 1 — max (cos (/)o, sin 0o) 



(15) 



(We treat K as a constant.) The solution of Eq.(jT^) with 
— / ^ a.t X — Q, and the cos and sin terms contribut- 
ing separately in the regions a; > 0, x < 0, is 



00 (a^) 



2 cos ^ tanh (u — a;/t«) , a; < 

|- — 2 cos^"' tanh (m + x/w) , a; > 



(16) 



where w = {K/A-kX'Y^'^ = hv/Ach, u = tanh ^(cos7r/8). 

The problem for the fiavor fields is simplified because 
they vary in the region where 0o ~ 7r/4. This gives 



|0'2 = 23/2^A'(1- 



cos 



sm 



(17) 



where cj) = 0i.2,3- Integrating Eq.(O) we obtain the func- 
tion 0(x). As illustrated in Fig.|l|7thc 0o step is ^ Vk 
times wider than the 0i,2,3 step. We emphasize that com- 
posite solitons in which all four fields vary are charge 
excitations of the lowest possible energy. 
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FIG. 1. Electron soliton charge and flavor parts 

[Eqs.(|l6|),(0),(|l|)] scaled by n/2, with the stiffness K = 16. 
The scales for (f)o{x) and (i>\.2:i{x) differ by K^^'^. The ap- 
proximate and exact solutions (|l9^, ( |l^ ) agree to 0.2%. 

Now we consider multi-soliton solutions and deter- 
mine charge compressibility from inter-soliton interac- 
tion. The compressibility x = i<PEn/(Pn)~^, with n 
the electron density, is directly related to the capacitance 
C = dn/dVg and the charging spectrum measured in a 
Coulomb blockade experiment. The charging spectrum 
peak spacing is 5„ = En+i - 2En + En-i « ■ We 
consider densitites smaller than Zj^ for which the charge 
fields 00 of different solitons may overlap, while the neu- 
tral cores with steps of 0i,2,3 are isolated. The effec- 
tive energy of the field 0o is obtained by minimizing (|^) 
with respect to 0i,2.3 at each value of 0Oi which gives 
— 4A' maxm cos ^0o + f ■ Switching between different 

branches occurs at 0o = f + f^^'- Approximating cos by 
a parabola near each maximum one obtains 



U{(j)o) = 2A' min„i ( 0o ^ nm/2 



(18) 



This approximation is extremely accurate, as illustrated 
in Fig.ljby comparing Eq.([l^) with the soliton solution 



x<Q 
X > 



(19) 



of the variational problem for the total energy 

E[^o\^ Jdx(^^dJoKdJo + U{4>o)^ (20) 

with the interaction kernel K = S{x — x') + ^V{x — x'). 

In a soliton lattice the field 0o {x) is a continuous mono- 
tonic function with smooth tt/2 steps, as in Figj^. To 
evaluate the energy (|2C|), we introduce a periodic dis- 
continuous field if — (fiQ — 7rm/2 with integer m cho- 
sen to minimize J7(0o)- The energy ( |20| ) with 0o = 
f ~ ^ J2m ^(^ ~ '^'*) takes the form 
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I dx 



^(a,(^+fG(^))i^(a,(^+fG(x))+2AV' (21 



where G{x) = ^j^S{x — ma). The advantage of the 
form (^) is that this quadratic function can be easily 
minimized in the Fourier representation. We obtain 

^(q) - -iq^G{q)K,/ {K,q^ + AnX') (22) 

with G{q) = 27r^jj S{qa — 27m). The energy of ( p^ ) is 

=5]4A%(fG(q))V(i^«9'+4^A') (23) 
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Using the identity (q) = L G{q) , with L the system 
size, the energy density can be written as 



E{n) = Tt.v—— 
8 



rn— — oo 



47rA' 



with qm = 2TTm/a and sohton density n = I /a. 



(24) 



electric field is screened and confined to ID at r < w and 
becomes deconfined at r > w. 



The crossover density n : 



is sensitive to the mag- 



netic field which controls soliton size w = hv/Ach- Thus 
the compressibility dependence on soliton size (Fig.^) is 
manifest in the charging spectrum dependence on mag- 
netic field. Along with the power law field dependence 
(^) of the gap at half-filling, it represents novel ID elec- 
tron correlation phenomenon observable in a narrow gap 
state of nanotubes in thermodynamic equilibrium. 
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2. Charge 
ity X ~ {d? E{n) / (P'n)^^ scaled by 1/hv plotted for several 
soliton widths w — (47rA')^^''^. Eq.(M) was used with the 
screening length Is — 0.2 /im, v — 9, ■ lO'^cm/s, e — 12, tube 
diameter d = 2r = 1 nm. 

Charge compressibility obtained from Eq.(^J) with 
2 4e2 



K{q) = 1 



1 TtflV 



\n[{q^ + T-^)/{q^+l-^)\ (25) 



that models screening (e.g., by a gate) at distance Ig ^ r 
is plotted in Fig.|^. At high density n 3> w^^ , compress- 
ibility is density- independent, Xo^ — f + tIt li^Cs/^)- 

At lower n <^ it varies as — Xo^ ~ tIt lii(?'ft). 
The crossover can be explained by noting that neighbor- 
ing overlapping solitons interact via V{r) cx r, while non- 
overlapping solitons interaction is V{r) oc 1/r, because 
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